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The family of q-Gaussian and q-exponential probability densities fit the statistical behavior of 
diverse complex self-similar non-equilibrium systems. These distributions, independently of the 
underlying dynamics, can rigorously be obtained by maximizing Tsallis “non-extensive” entropy 
under appropriate constraints, as well as from superstatistical models. In this paper we provide an 
alternative and complementary scheme for deriving these objects. We show that q-Gaussian and 
(/-exponential random variables can always be expressed as function of two statistically independent 
Gamma random variables with the same scale parameter. Their shape index determine the com¬ 
plexity g-parameter. This result also allows to define an extended family of asymmetric g-Gaussian 
and modified g-exponential densities, which reduce to the previous ones when the shape parameters 
are the same. Furthermore, we demonstrate that simple change of variables always allow to relate 
any of these distributions with a Beta stochastic variable. The extended distributions are applied 
in the statistical description of different complex dynamics such as log-return signals in financial 
markets and motion of point defects in fluid flows. 

PACS numbers: 02.50.-r, 89.75.Da, 89.65.Gh, 47.27.-i 


I. INTRODUCTION 


Long-range interparticle interaction, long-term micro¬ 
scopic or mesoscopy memory, fractal or multifractal oc¬ 
cupation in phase space, cascade transfer of energy or 
information, and intrinsic fluctuations of some dynam¬ 
ical system parameters are some of the properties that 
nowadays are related with complexity. One of the emer¬ 
gent properties related with these phenomena is the 
power-law statistics of the corresponding nonequilibrium 
states. While there exits different underlying formal¬ 
ism for tackling these issues, maximization of Tsallis 
“non-extensive” entropy provides an alternative ba¬ 
sis over which complexity can be analyzed and studied 
in a broad class of systems. Introducing a generalized 
second moment constraint ii , the formalism lead to 
a generalization of standard normal probability densi¬ 
ties, named as g-Gaussian distributions. In terms of a 
generalized exponential function, they can be written 
as P{x) = {\fP/Mq)eyip>q{—()x^), where the parameter 
g G Re defines different complexity classes. Explicitly, 
these statistical objects read 


P(x) 


2? 


I - (I - g)/3x^ 



—00 < g < I, (I) 


where the variable of interest x is restricted to the domain 
0 < (1 — q)fix^ < I. On the other hand, 


P{x) = 


VI 

Mq Ll -b (g - l)j)x‘^. 


1<9<3, (2) 


of the distributions. As is well known Q, in the limit 
g —>■ I both expressions reduce to the standard Gaus¬ 
sian distribution. These generalizations allow to describe 
variables restricted to a finite domain [Eq. ([T])] as well as 
power-law statistics [Eq. (2))]. 

g-Gaussian distributions also arise as solution of non¬ 
linear Focker-Planck equations [1, @ as well as in the 
formulation of central limit theorems with highly corre¬ 
lated random variables Q- Furthermore they fulfill a 
generalized fluctuation relation symmetry Q. A wide 
and diverse class of systems obey their statistics Q , such 
as in fluids flows [9l-[lll| , optical lattices [l^ , trapped ions 
interacting with a classical gas , in granular mixtures 
anomalous diffusion in dusty plasma m or cellu¬ 
lar aggregates [l^, avalanches sizes in earthquakes m, 
in astrophysical variables [l^ , as well as in econophysics 
|l9l - i2^ . 

When introducing a first moment constraint, Tsallis 
entropy leads to a g-exponential distribution [Eqs. ([1]) 
and m under the replacement —>■ x with — oo < g < I 
and I < g < 2 respectively], which in the limit g —>■ I re¬ 
covers the standard exponential probability density of a 
positive random variable. These objects, for example, al¬ 
low to fit high energy collisions [2^ , quark matter statis¬ 
tics [ 2 ^, solar flare s and momentum distributions 
of charged hadrons |27j| . g-exponential functions also fit 
anomalous power-law dipolar relaxation as well as 
spin-glass relaxation [ 2 ^ . More recently, a kind of gener¬ 
alized g-Gamma probability density was introduced for 
describing stock trading volume flow in financial markets 


where now x is allowed to runs over the real line. The re¬ 
striction g < 3 follows from the normalization condition 
dxP{x) = I, which is guaranteed by the dimension¬ 
less constant Afq. The parameter y/P measures the width 


It is remarkable that all quoted probability densities 
can also be obtained from a superstatistical modeling Q , 
where a parameter of an underlying probability measure 
becomes a (positive) random variable characterized by 
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a Gamma distribution [ssl - l^ . This is the case of q- 
Gaussian densities, where the underlying distribution is 
a normal one [3, while for g-exponential it is an expo¬ 
nential function [^. For generalized g-Gamma variables 
the underlying distribution is a Gamma density while the 
random parameter is distributed according to an inverse 
Gamma distribution [s^. 

The main goal of this paper is to present a complemen¬ 
tary and alternative scheme to those provided by entropy 
extremization and superstatistics. We show that random 
variables described by any of the quoted families of g- 
distributions can be written as a function of two indepen¬ 
dent (positive) random Gamma variables Their 

scale parameter is assumed the same, while their shape 
indexes determine the complexity g-parameter. When 
the shape indexes are different, a class of extended asym¬ 
metric g-Gaussian and modified g-exponential distribu¬ 
tions are obtained. Generation of g-distributed random 
numbers is straightforward from these results (s^ . [s^. 
We also show that simple transformation of variables al¬ 
low relating any of the obtained densities with a Beta 
distribution. Interestingly, this statistical function has 
been applied to model a wide variety of problems arising 
in different disciplines On the other hand, a g- 

triplet [ 3 ^, 113 for the probabilities densities is obtained. 
The usefulness of the extended distributions in the con¬ 
text of financial signals [2l| and motion of point defects 
in fluid flows [TTj | is demonstrated. These systems are 
characterized by highly asymmetric distributions. This 
feature is absent in previous approaches, being recovered 
by the present one. 

The paper is organized as follows. In Sec. II we review 
the properties of Gamma random variables and introduce 
the main assumption over which the present scheme re¬ 
lies. In Sec. Ill asymmetric g-Gaussian distributions 
are obtained, while Sec. IV is devoted to modified g- 
exponential densities. In Sec. V the properties of the 
proposed scheme as well as its applications are discussed. 
In Sec. VI the Conclusions are provided. 


II. MODEL 

A stochastic random variable Y is Gamma distributed 
if its probability density is 

, e-y/^ 

P{y)=y°‘ 0aY{a) ' 0<a<oo, (3) 

where r(a) is the Gamma function. This distribution 
is characterized by the scale parameter 9 > 0 and its 
shape parameter a. In terms of these parameters its 
mean value reads (Y) = yPiy)dy = a9, with variance 
var[y] = (Y^) — (V)^ = a9'^. The underlying stochastic 
process that leads to this statistic involves a cascade-like 
mechanism In fact, this property is evident from 

the Laplace transform P{u) = -I- where 

P{u) = dtP{t)e~'^*. Hence, when a is natural it re¬ 


duces to a convolution of exponential functions, which 
can be read as a cascade of consecutive random steps. 

The present approach relies on two independent 
Gamma random variables Yi, ^ 2 - Their joint probabil¬ 
ity density then read 


Pi3P.y2) = yr^ 



g-(yi+V2)/6 

0“+“T(a)r(a')' 


( 4 ) 


Here, we assumed that both scale parameters 9 are the 
same, while a and a' are the shape parameters of Yi 
and Y 2 respectively. The main ingredient of the present 
scheme is the ansatz 


X = f{Y„Y2), (5) 

where the new random variable A, depending of the func¬ 
tion f{yi,y 2 ), develops different statistics. We will show 
that a wide class of g-distributions arise from non-linear 
functions, which in turn are asymmetric in their argu¬ 
ments (see Sec. HI and IV). Nevertheless, in all cases 
they fulfill a very simple symmetry (see Sec. V). 

The probability distribution of X is completely de¬ 
termined by the joint probability (U) and the function 
f{yi,y 2 )- In fact, it follows from a elementary change of 
variables [3^. For closing the problem, we introduce an 
extra random variable Z defined by the addition 


Z={Y,+Y2). (6) 


Therefore, the joint probability of X and Z is given by 


P{x,z) = P(j/i,j/ 2 )|det(J)|, 
where J is the Jacobian matrix 


J = 


dyi 

dyi 

dx 

dz 

dyi 

dyi 

dx 

dz 


( 7 ) 


( 8 ) 


The probability of each variable follows by partial inte¬ 
gration 


P{x) 



dzP{x, z), 


p{z) = r 

J —c 


dxP{a 


(9) 


As Z is defined by the addition of two independent 
Gamma variables with the same scale factor, it follows 


P{z) 


_ —1 


g-ve 

6»“+“T(a + a')' 


( 10 ) 


Hence, Z is also a Gamma variabl^Z > 0, see Eq. ©] 
where its shape index is {a + a') 


III. ASYMMETRIC g-GAUSSIAN 
DISTRIBUTIONS 

In order to motivate the election of the function 
/(yi: J/ 2 ) that lead to g-Gaussian statistics we may think 
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(in a roughly way) in a Brownian particle that inter¬ 
act with a complex reservoir. Yi and Y 2 are the (posi¬ 
tive and negative) impulse moments induced by the bath 
fluctuations. In addition, the complexity of the system- 
environment interaction is taken into account by a sys¬ 
tem response function M“^(Yi, 1 ^ 2 ) that depends on both 
Yi and Y 2 . Therefore, this contribution can be read as 
a random-dissipative-like mechanism. The particle fluc¬ 
tuation is finally written as X Ri (Yi — Y 2 )/M{Yi,Y 2 ). 
One may also think in a economical agent that, from the 
available information, predict that a future price may in¬ 
creases or decrease a quantity Yi or Y 2 respectively. The 
weight of the available information is then measured by 
M“^(Yi, Y 2 ), leading to the same kind of dependence. 

In order to close the model, we assume that M(Yi, Y 2 ) 
is given by a kind of average or mean value between the 
two random values Yi and Y 2 . Specifically we take 


M(Yi,Y2) 


lirr + y^) 


( 11 ) 


where /i G Re. Therefore, we write the X random variable 
[Eq. ([1])] as 


X = 


El - Y 2 




^{Yr+Y^^) 


1/m ' 


( 12 ) 


By convenience we introduced a factor one half. On the 
other hand, the additional parameter > 0 scales and 
gives the right units of X. In fact, notice that the re¬ 
maining contribution in Eq. © is dimensionless. Tak¬ 
ing different values of the real parameter fi a wide class 
of probability distributions arise, which in turn may also 
depends on the parameters 0, a and a' that determine 
the joint probability density ([4]). 


A. Arithmetic mean 


The arithmetic mean value corresponds to /i = 1, im¬ 
plying that 


1 Y 1 -Y 2 
y?Ti+Y2' 


(13) 


Notice that, for any possible value of Yi and Y 2 , the ran¬ 
dom variable X assume bounded values in the domain 
{—l/y/P,+l/y/P). Taking into account the Z variable 
[Eq. ([0])], we obtain the following inverted relations 


Yi = |(l + v^X), Y 2 = ^{l-./pX), (14) 


which in turn implies that | det(J)| = Eqs. dH) 

and © lead to P{x, z) = P{x)P{z), where P{z) is given 
by Eq. Therefore, the random variables X and Z 

are statistically independent. Eurthermore, X obeys the 
statistics given by the probability density 

p{x) = ^{i + ^xr-\i-^xr'-\ ( 15 ) 

OlOl' 


where the normalization constant reads Maa' = 
2“+“'“ir(Q!)r(a')/r(a -I- a'). Notice that P{x) does not 
depends on the scale parameter 9 [see Eq. (jd])]. It only 
depends on the shape indexes a, a', and the scale param¬ 
eter /3. 

The distribution ()15p develop a maximum located at 


xm = 


a — a 


V? (a + a' - 2) ’ 


(16) 


when a > 1, a' > 1, or at xm = ±l/v(d in any other 
case. Its average value read 


(X) = 


1 a — a' 


y/P aY a'' 

while the variance var[X] = (X^) — (X)^ is given by 


(17) 


var[X] = — 


daa' 


P {a + a')^(l -bad- a') 


(18) 


q-Gaussian distributions 


Eq. (fT^ can be rewritten as 


Hence, we name this function as an asymmetric Poisso- 
nian g-Gaussian distribution G^^{x\q, a,/3) with index g, 
and asymmetry parameter a. 




( 20 ) 


Erom the positivity of a and a' the asymmetry index 
must to satisfy |a| < (2 — g)/(l — g). In the symmetric 
case, a = 0, a' = a, we get 


P(cr) = ^(l-^x2)“-\ (21) 

where Afa = 2^““^r^(a)/r(2a). Therefore, under the as¬ 
sociation /3 —>■ /3(1 — g), with a > 1, we recover Eq. ([T]). 
Over the domain a G (l,oo), the non-extensive parame¬ 
ter runs in the interval q G (— 00 ,1). 

In Fig. 1 we plot the function (I19p for different values 
of the asymmetric factor a, Eq. ([201). For increasing 
a > 0, the distribution accumulates around y/Px ~ 1. 
For a < 0 a reflected accumulation around y/px r; — 1 is 
developed. The distribution with a = 0 corresponds to 
Tsallis non-extensive thermodynamics. 


B. Geometric mean 


In Eq. dill) the geometric mean correspond to lim/r —>■ 

r 1 1//J. 

0, which satisfies lim^^o 5 ( 1 ^ + ^ 2 ) — 
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FIG. 1: Poissonian g-Gaussian probability distribution Eq. 
HU}, for different values of the asymmetry a, Eq. pO]) . The 
circles correspond to a numerical simulation based on Eq. 

m- 





FIG. 2: Poissonian g-Gaussian probability distribution Eq. 
(|30|l . for different values of the asymmetry a, Eq. (I31|l . The 
circles correspond to a numerical simulation based on Eq. 

dig). 


Therefore, we get the random variable [Eq. (11^) ] 


1 Y 1 -Y 2 


( 22 ) 


Notice that X assumes values over the entire real number 
line, X S Re. In this case, the inverted relations are 


Due to the previous asymptotic behaviors the moments 
are not defined for any value of the characteristic shape 
parameters. When a > 1/2 and a' > 1/2, the average 
value reads 




„ r(a-i)r(a'-i) 
“ ’ 2 r(a)r(a') 


(28) 


El = 


f( 


1 + 


^X 


+ /3E2 


), U.f(l- 


^X 


(23) 


implying that |det(J)| = y/P{z/2){l -\- jix^) which 
in turn also lead to P(x, z) = P{x)P{z), but here 


P{x) = 




/ V1 + Bx'^ ) 


Afaa' V 1 + Bx'^ 


(■\/l + Bx"^ + ^/px) 


(\/i + Bx'^ - s/BxT 


(24) 


As in the previous case, this distribution is independent 
of the rate parameter 9. The normalization constant is 
the same, Maa' = 2“+“ “^r(a)r(a')/r(a + a'). 

Eq. (1241) develops a maximum, which occurs at 


XM 


a — a' 

\/(l + 2 q!)(1 + 2q,'')B 


(25) 


In the limit \/Bx ^ 1, a power-law behavior arise 


lim P{x) 


PP na-c 


AfOta' 


\y/BxJ 


2a'+ 1 


(26) 


while for y/Bx <C — 1 we obtain 


while the second moment, for a > 1 and > 1, is 


)=fl[(a-a) +[a+a-2)\ - , . , -. (29) 

B 41 (a)l (a') 


Outside the previous intervals the first two moments are 
not defined. 


q-Gaussian distributions 


Eq. (l24ll can be rewritten as 


P{x) 



{ ^ 1 

a + a' + l 

1 ^ 1 

/ -\/i -1- Bx^ + VBx\ 

•^OLO.' 

V 1 -1- 

^ -\/i + Bx"^ — VBx' 


(30) 

As is the previous case, we name this function as an asym¬ 
metric Poissonian g-Gaussian distribution G^i{x\q,a, B) 
with index q, and asymmetry parameter a. 


g=l + 


'a + a' 


L 2 




(31) 


Hence, 1 < g < 3 and the asymmetry factor satisfy the 
restriction jaj < (3 —g)/[2(g— 1)]. In the symmetric case, 
a = 0, a' = a, P{x) reduces to 


lim P{x) 

x^ — oo 


PP aa’-a( 1 A 
Maa' \-y/Bx) 


(27) 


P{x) 


VP ( 1 

Afa ^ 1 + Px'^ ) 


(32) 









































5 


where Me = 2^““^r^(a)/r(2a). Under the association 
/3 —>■ f3{q — 1), we recover Eq. ([5]). In the interval a S 
(0, oo) the non-extensive parameter runs in the interval 
qe (3,1). 

In Fig. 2 we plot the function (l30l) for different values 
of the asymmetric factor a, Eq. dSlD. For increasing a > 
0, the maximum of the distribution is shifted to higher 
values. For a < 0 the extremum develops for negative 
values. The symmetric case a = 0 corresponds to the 
q-Gaussian distribution arising from Tsallis entropy. 


Notice that this assumption is very similar to Eq. m, 
but here X is a bounded positive stochastic variable, 0 < 
y/PX < 1. Using the approach defined in Sec. II, here 
we obtain the inverse relations 

Yi=z{l-y/px), Y2=Zy/pX, (36) 

while I det(J)| = y/^z, leading again to a statistical inde¬ 
pendence of X and Z, that is P{x,z) = P{x)P{z). The 
density of interest here is 


C. Relation between both cases 


P{x) = -^i^xY 

Vrva' 


- 1(1 - y/YxY 


(37) 


Given X determinate by relation (II3, the random 
variable X' defined as 


X' 


X 

x/i - YX^ ’ 


(33) 


recover Eq. (HH). This simple relation demonstrate that 
there exist a one to one mapping between g-Gaussian 
variables in the different domains of the complexity pa¬ 
rameter q. In fact, if we defmBy/pX = sin((/)) € (-1,1) 
for q G (— 00 , 1 ), hence yfpX' = tan(((i) G (—oo,-|-oo), 
where X' has associated the index q G (1,3). 

Alternatively, if X is given by Eq. (12^ . the inverse 
transformation 


X' 


X 

v'l + px^ ’ 


(34) 


lead to Eq. m- While these relations are known for 
symmetric g-Gaussian distributions (a = 0) Q, here we 
showed that they are also valid for the asymmetric den¬ 
sities (a Y 0) introduced previously. 


IV. MODIFIED g-EXPONENTIAL 
DISTRIBUTIONS 

Variables distributed according to a g-exponential den¬ 
sity are positive. Therefore, the previous scheme does not 
apply, but a similar one can be implemented. We name 
the emerging distributions as modified g-exponential den¬ 
sities. Taking into account the notation of Refs, (sol - ls^ 
they can also be called as generalized g-Gamma densities. 
Nevertheless, it seems that they do not satisfy the same 
entropic properties than standard Gamma distributions 
(^ . On the other hand, we remark that some proper¬ 
ties of the following distributions are known and can be 
found under different denominations [sj, ■ 


A. Bounded domain 


For getting a positive variable, we introduce the fol¬ 
lowing functional dependence 


X = 


1 U2 

yYPY^+Y2' 


(35) 


where Maa' = r(a)r(a')/r(a -|- a')- When a > I and 
a' > I, P{x) reaches a maximal value located at 


I a' - I 
^(a + a'-2)^ 


Its first moment reads 


(^) 


1 a' 

y/p aP a'' 


while the variance is given by 


var[X] 


1 aa' 

P {a + a')^(I P a + a') 


(38) 


(39) 


(40) 


q-exponential densities 

The distribution (EZl) may be named as a modified 
Poissonian g-exponential distribution E^]^(a;|g, d, P) with 
index g, and “distortion parameter” d, 

g = 1-, d = a' — 1. (41) 

a — I 

Therefore, —oo < g < I and d > —I. When d = 0, that 
is a' = I, Eq. dSll) becomes 

P{a:) = ^{l-y/Pxr-\ (42) 

where Afa = r(a)/r(a -b 1). Therefore, under the ex¬ 
tra association y/p —>■ y/Pil — g), we obtain a stan¬ 
dard g-exponential density. For a G (I,oo), it follows 
g G (- 00 ,1). 

In Fig. 3 we plotted the function dazi) for different 
values of the distortion parameter d, Eq. SB). For d < 0, 
the density diverges around the origin. This property is 
inherited from the Gamma distribution Eq. (|3]). On the 
other hand, for d > 0 the density vanishes at the origin 
and for increasing d it accumulates around y/px ~ I. The 
plot for a = 0 is the g-exponential distribution arising 
from Tsallis entropy. 
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FIG. 3: Poissonian g-exponential probability distribution Eq. 
(|33, for different values of the distortion parameter d, Eq. 
(HU- The circles correspond to a numerical simulation based 
on Eq. (l35l) . 



FIG. 4: Poissonian g-exponential probability distribution Eq. 
(|45|l . for different values of the distortion parameter d, Eq. 
poll. The circles correspond to a numerical simulation based 
on Eq. (l43ll . 


B. Unbounded domain 


while the second moment is only defined for a > 2, 


An unbounded positive variable (0 < AT < oo) is ob¬ 
tained from the relation 




(43) 


ir(a- 2 )r(a'-b 2 ) 
^ r(a)r(a') 


(49) 


where as in the previous cases scales the random q-exponential densities 

variable X. Here, the inverse relations are 


Yi = 


Z 

1 -i- ^/px 


^ z^x 

" l + ^/px' 


(44) 


while |det(J)| = -I- , leading to P(x,z) = 

P{x)P{z), where 


P{x) 


yfP iVPx)Y 1 

A/'aa' (1 + y?a;)“+“' 


(45) 


with J\faa' = r(a;)r(a')/r(a -I- a'). This distribution de¬ 
velops a maximum that is located at (a' >1) 


The distribution (1371) may also be named as a modified 
Poissonian g-exponential distribution £’^]^(x|g, d, /3) with 
index q, and distortion parameter d, 

q = 1 -\ - d = a' — 1. (50) 

a -|- q; 

In consequence, 1 < q < oo, and the distortion parameter 
satisfy (2 — g)/(g — 1) > d> —1. When the distortion is 
null, a' = 1, Eq. (|45l) becomes 


Xm 


1 a'-1 

\/?(a + l)’ 


(46) 


P(x) 


VP 1 
A4 (1 -b \/?a;)“+i ’ 


(51) 


FoiVP^ » 1, it follows the asymptotic power-law behav- 
ior 


lim F(x) 



(47) 


In consequence, the moments are not defined for any 
value of the shape parameter a. For a > 1, the mean 
value reads 


1 r(a-l)r(a' + l) 

^ ^ ^ T{a)r{a') 


(48) 


where Afa = r(a)/r(l-|-a). Thus, under the extra associ¬ 
ation y/P —>■ VPV~q) we get a g-exponential probability 
density. In this case, for a S (0, oo), it follows g G (2,1). 

The function (j45|) is plotted in Fig. 4 for different 
values of the distortion parameter d, Eq. (IhHI) . For d < 0, 
the density diverges around the origin. On the other 
hand, for d > 0 the density vanishes at the origin. In all 
cases a power-law behavior is obtained for y/px » 1. The 
plot for d = 0 is the g-exponential distribution arising 
from Tsallis entropy. 
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C. Relation between both cases 


Given X determinate by Eq. (1351) . the random variable 
X' defined as 


X' = 


1 - 


(52) 


is given by Eq. (03). Alternatively, if X is given by Eq. 
(03- the inverse transformation 


A' = 


A 

l + y?A’ 


(53) 


lead to Eq. (1551) . These conjugate relations are valid 
for both the unmodified {d = 0) as well as the modified 
{d 7 ^ 0) ^-exponential densities. 


D. Stretched g-exponential densities 


Levy distributed numbers was also established [s^. On 
the other hand, numerical generation of Gamma random 
numbers is also well established [13, • Therefore, the 

present scheme defines an alternative and solid basis for 
obtaining ^-distributed random numbers by generating 
two independent Gamma random numbers. Using this 
method, in Fig. (1) to (4) we explicitly show (circles) 
the recovering of the symmetric and unmodified distri¬ 
butions, all of then corresponding to Tsallis entropy for¬ 
malism. 


B. Symmetries 

While the underlying joint statistics of the Gamma 
variables depends on the scale parameter 9, Eq. the 
distributions of X do not depend on it. This is not the 
only symmetry of the proposed scheme. In fact, it is 
simple to check that the symmetry 


Introducing the change of variables 

^Jjx=[^x)^'\ (54) 

defined by the extra parameter v G Re, if X is distributed 
according to Eq. ([37|) . it follows [P{x)dx = P{x)dx] 



On other hand, if X is distributed according to Eq. (1^ , 
we get 


P{x) = 



(56) 


In both cases, imposing the condition va' = 1, the previ¬ 
ous two expressions becomes stretched g-exponential den¬ 
sities, P{x) = (-v/?/A/g) expg[-(v^x)'"] {x > 0, ly > 0). 
Hence, these distributions can also be covered with the 
present approach [Eqs. (1551) and (H51) under the change 
of variables (15^ ]. 


V. PROPERTIES AND APPLICATIONS 

In the previous two sections we demonstrated that the 
assumption allow us to recover and to define an ex¬ 
tend family of g-Gaussian and g-exponential densities. 
Here, we discuss some general properties of the approach 
as well as some applications of the extended distributions. 


A. g-distributed random numbers 

Numerical generation of random numbers obeying g- 
Gaussian statistics was explored previously by introduc¬ 
ing a generalized Box-Muller method [s^l ■ Generation of 


/(ih,r2) = /(^,^), (57) 


is fulfilled, where /(Ti,l 2 ) define the X random vari¬ 
able, Eq. ([S])- In fact, this property is valid for the 
g-Gaussian case [Eq. (IT^ for any /r] as well as for the 
g-exponential variables [Eqs. (1^ and (l4^ ]. We notice 
that f{Yi,Y 2 ) = g{Yi/Y 2 ), for arbitrary functions g{y), 
always satisfies the relation (1571) . Extra structures can be 
established by introducing arbitrary change of variables. 

The symmetry (|57l) implies that the same results arise 
if instead of Gamma distributed variables one take in¬ 
verse Gamma variables, that is, Y' = 1/Y where Y is 
Gamma distributed [Eq. ([3])]- Using that P{y)dy = 
P{y')dy', it follows 


Piy') 


1 

(y/)a-n go^Y{a)’ 


0 < a < oo. (58) 


C. Relation with Beta distributions 

As all functions f{Yi,Y 2 ) fulfill the condition (1571) . it 
is clear that any of the corresponding variables X are 
always related by a change of variables between them. 
Therefore, it does not make sense to affirm that one of 
them generates or is more fundamental than the others. 
Nevertheless, here we want to emphasize that any of the 
probability densities obtained previously can be related 
with the well known Beta distribution (33l - [3^ . It reads 

I (q;)I [a ) 

where the domain of its variable is w G (0,1). Further¬ 
more, its shape parameters a and a' are positive. 

Defining the change of variables w = w{x), all obtained 
g-distributions becomes equal to Eq. (1551) . Alternatively, 












defining a new variable x = x(w), from the Beta distri¬ 
bution it is possible to obtain the g-densities. Explicitly, 
for the asymmetric distribution Eq. (ITKI) [or Eq. (IISl)], 
the change of variables read 

= ^(1 + ^x = 2(w-^). (60) 

Therefore, all (asymmetric and symmetric) g-Gaussian 
distribution with — oo < g < 1 are related to a Beta 
variable by a shifting of their arguments. 

For the g-Gaussian defined by Eq. (l24l) [or Eq. W], 
where 1 < g < 3, the change of variables is 


any of the extended distributions can be rewritten as 
P(x) = (v^/A/’q) exp^[—/?E(a;)]. This structure emerges 
from Tsallis entropy by using a constraint based on a gen¬ 
eralized mean value of V(x) Q. Nevertheless, here the 
resulting functions V (x) depend on the parameter g and 
also on the asymmetry and distortion factors. Therefore, 
a relation between Tsallis entropy and the asymmetric 
and modified distributions cannot be established in this 
way. 

D. q-triplet for probability densities 




\/?x \ 

-I- /3x^'' ’ 


= 



x/w(l - w) 


• (61) 


For Eq. (EZD, the relations are 


w = 1 — \/Px, = 1 — w, (62) 


that is, the modified (and standard) g-exponential den¬ 
sities in the interval — oo < g < 1 arise from an axe 
inversion of a Beta distribution. Finally, in the interval 
1 < g < oo, Eq. (|151) . the transformations are 


w = 


1 -I- v(dx’ 


\/Px = 


I — w 


(63) 


The previous relations can be enlighten by using that 
a variable W obeying the Beta statistics (l5^ can also be 
written in terms of two independent Gamma variables 
(El and Fa) [Eq. ®] 


W = 


Vi 

El + Es ’ 


W' = 


V2 

El + E2 ’ 


(64) 


where the additional variable W' is also Beta distributed. 
In fact, W -I- W = 1. After a straightforward manipula¬ 
tion, the random variables associated to the g-Gaussian 
distributions, Eqs. m and (l2^ . can respectively be 
rewritten as 


= {W-W'), 




(65) 


while for the g-exponentials densities, Eqs. (l35|) and (l4^ . 
respectively it follows 


= W', (66) 

Both Eq. (1651) and Eq. (IMl) show the stretched rela¬ 
tion between all the generalized g-distributions and Beta 
random variables. In fact, any stochastic variable defined 
by a function satisfying the symmetry (IS3 can be related 
by a transformation of variables with a Beta distribution, 
Eq. dSH). 

One may also take the inverse point of view and to 
explore if the previous densities can be obtained from 
Tsallis entropy under a more general constraint. In fact. 


In the context of non-extensive thermodynamics three 
different values of the complexity parameter, named as 
g-triplet, are associated to different physical properties 
such as the statistics of metastable or quasi-stationary 
states, sensitivity to initial conditions, and time-decay 
of observable correlations [s^ H^. Here, we show that 
three different values of g allow to indexing the symmetric 
and unmodified probability densities. We remark that 
not any direct relation can be postulated between both 
triplets, because here it is established for normalizable 
objects, P(x)dx = 1. 

We denote by g® and g® the complexity indexes of 
the g-Gaussian distributions Eqs. (EOI) and m respec¬ 
tively. Furthermore, g^j and denote the indexes of 
the g-exponentials, Eqs. m and (1501) respectively. The 
four indexes are given by 


y<i 


9<i 


1 

1 


1 

a- 1’ 
1 

a — 1 ’ 


= 1 + 
9>i = 1 + 


1 


a-f 1/2 

1 


ct -)- 1 


(67a) 

(67b) 


We notice that g®^ = g®^ [see Eqs. (1001) and (HTI) ]. 
This equality is expectable because g-Gaussian and g- 
exponential distributions for — oo < g < 1 are related 
by a linear change of variables [see Eqs. dOO]) and (EH)] 
with a Beta distribution. Therefore, the complete family 
of analyzed distributions can be indexed with only three 
g-parameters: (g^j^, g® g^^^). The previous expressions 
are equivalent to 


1 

1 

l-9<i 


q; — 1, 
a — 1, 


1 1 
g®i-l “”^2’ 


9>i-l 


= a + l. 


(68a) 

(68b) 


From here we realize that there exist simple relations 
between any of the g-triplet parameters. 


E. Applications of the extended distributions 

The assumption ([5]) lead us with a broad class of prob¬ 
ability densities, which in turn cover the most used prob¬ 
abilities densities arising from Tsallis entropy maximiza¬ 
tion. Hence, besides it theoretical interest, we ask about 
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the possible applications of the asymmetric and modified 
distributions. 

From the previous analysis, we arrived to the con¬ 
clusion that asymmetric g-Gaussian and modified q- 
exponential distributions in the interval —oo < g < 1 
[Eqs. (IT^ and ([37)l ] are related by a linear change of vari¬ 
ables with a Beta distribution. Therefore, these functions 
fall in the wide range of applicability of this distribution 
[33l - l3^ . For example, (unnormalized) Beta distributions 
emerge in the statistical description of quark matter (see 
Eq. (102) in Ref. El). 

The modified g-exponential function Eq. (HSl) was used 
in the description of stock trading volume flow in fi¬ 
nancial markets E&ii (named as generalized g-Gamma 
probability density). This distribution is also known as a 
Pearson Type VI distribution or alternatively Beta-prime 
distribution El (see also El)- 

To our knowledge, asymmetric g-Gaussian distribu¬ 
tions Eq. (1501) were not used previously. In the present 
approach, the asymmetry of these probability densities 
has a clear dynamical origin. In fact, associating a cas¬ 
cade process to each Gamma variable, asymmetries arise 
whenever the cascades have a different shape index (a 
and a'). Below we discuss the application of these kind 
of distributions as a fitting tool in the context of financial 
signals El ^'iid movement of defects in fluid flows El. 


1. Log-returns signals on large time windows 


Erom the price signal y{t) in a financial market it is 
possible to define the stochastic process y{t) = ln[?/(t -1- 
At)/y{t)], where At is a constant time interval. This 
signal gives a simple way of representing returns in the 
market. Usually it is studied the normalized log-returns 
ZAt{t) = ln[y(t) - {y{t))]/<JAt, where {y{t)) is the av¬ 
erage and aAt gives the standard derivation of y{t) for 
a given At. Daily closing price values of the S&P index 
for a period of twen ty y ears were analyzed by Ausloos 
and Ivanova in Ref. |2l| . Assuming a stationary signal, 
ZAt(t) —)> ZAt, for large time windows {At > 1 day), the 
authors fitted the experimental data with a g-Gaussian 
like distribution 



where {1/Afg) = Q;r(^)/[r(^ - ^)F(^)] and 
depends on the parameters q and a (see Eqs. (4) and (5) 
in El). This distribution can be obtained from a super- 
statistical model assuming, for example, that Brownian 
particles diffuse in a potential U{x) = C\x\'^°- On 
the other hand, Eq. (l6^ can also be recovered from the 
present approach based on random Poisson variables. In 
fact, it follows by extending symmetrically (x —>■ |a:|) the 
stretched g-exponential distribution (I56p , with a' = \/v, 
and the following replacements x —>■ \zAt\i $ ^ fdq, 
V —>■ 2d, and a —>■ (^^ — ^). Random numbers gen¬ 
eration is achieved by introducing an extra stochastic 




At (days) 



At (days) 

FIG. 5: (a) Probability density p{zAt) of normalized log- 
returns for the S&P index (circles). The full line corresponds 
to Eq. (1241) (see text). The inset shows the peak region. The 
parameters are a = 1.97, a' = 2.85, and /3 = 0.37. From (1311) 
it follows q = 1.34 and a = —0.87. (b) Dependence with At 
of the shape parameters a and a'. Dotted line linear fit. (c) 
Parameter /3. 


variable that with probability one half defines their sign 
(positive or negative). 

Eq. ((551) develops an asymptotic (y^l^Atl 1) 
power-law behavior. Nevertheless, the authors also find 
that the experimental data are not consistent with the 
symmetry p{zAt) = p{—ZAt)- In particular, the power- 
law behaviors have different exponents for positive and 
negative values. Similar asymmetries were found previ¬ 
ously in Ref. [d^ . 

Here, we may associate the observed asymmetry of the 
data to two cascades mechanisms, each one being repre¬ 
sented by a Gamma random variable. In a roughly way, 
the difference between both variables can be associated to 
















10 


different networks properties related to the propagation 
of information that support an increased or decreased fu¬ 
ture value. The complete system response is defined by 
Eq. (22), that is a geometric mean value of the driving 
fluctuations. Hence, instead of using Eq. (l6^ . we pro¬ 
pose to fit the probability density of the log-returns with 
the asymmetric g-Gaussian distribution (1301) under the 
shifting P{x) —>■ P{zAt + (X)), where (X) is given by 
Eq. (IMl). 

In order to check this proposal, here we analyze the 
daily closing price values of the S&P index for the 
period between Jan. 3, 1950 and Dec. 3, 2014, which 
provides a 16336 data base. In Fig. 5(a) we show the “ex¬ 
perimental” probability distribution (circles) for At = 30 
days. The data are clearly asymmetric, which in fact 
are fitted by Eq. (l30ll (full line). Its characteristic pa¬ 
rameters a, a', and /3 were determinate by minimizing 
the global error. Based on the quadratic global error 
J2iPexp-Ptheory)'^/Pexp, we checked that for a wide range 
of At the asymmetric distribution provides a better fit¬ 
ting than Eq. (162- 

In Fig. 5(b) we plot a and a' as function of At. In 
the limit At —>■ 0 the asymmetry vanishes (a ~ a'). Fur¬ 
thermore, in the plotted interval, both shape parame¬ 
ters increase linearly with At. For higher values of At an 
irregular-logarithmic-like grow behavior is observed (not 
shown). For At > 500 the distributions approaches nor¬ 
mal Gaussian ones. This limit is consistent with the grow 
of a and a' [see Eq. (I5T]) ]. On the other hand, we find 
that 1//3 also increases linearly with At, Fig. 5(c). This 
(diffusive) behavior is also found for intervals At minor 
than a day [^ . 

Extra analysis and ingredients are necessary for ex¬ 
plaining the linear behaviors shown in Figs. 5(b) and (c). 
On the other hand. Fig. 5(a) shows that the proposed 
probability density provides a reasonable and alternative 
fitting to that based on Eq. (E2, which in turn is able 
to capture the observed asymmetries. 


2. Defect velocities in inclined layer convection 

In Ref. m Daniels, Beck, and Bodenshatz studied the 
motion of point defects in thermal convection patterns in 
an inclined fluid layer (heated from below and cooled 
from above), a variant of Rayleigh-Benard convection. 
Due to the inclination the system is anisotropic. The (ex¬ 
perimental) probability distribution of the (positive and 
negative) defect velocities is different in the transverse 
(i, across rolls) and longitudinal (y, along rolls, uphill- 
downhill) directions. In the transverse direction the ve¬ 
locity {vx) can be fit with a symmetric q-Gaussian distri¬ 
bution {q ~ 1.4). Nevertheless, in the longitudinal direc¬ 
tion (vx) the distribution, depending on a dimensionless 
temperature e (see details in [il|), develops strong asym¬ 
metries. In that situation, Tsallis distributions, even de¬ 
fined with a cubic potential, are unable to fit the experi¬ 
mental data (see Fig. 4(a) in [HI)- As shown in the next 


10° 
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CL 
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y 

FIG. 6: Probability density p(yy) of the normalized velocity 
Py = Vy!^ (aj) — {vy)'^ for positive and negative defects in 
inclined layer convection, e = 0.08 (see m for details). The 
fitting (dotted and dashed lines) correspond to Eq. (1301) (see 
text). For positive defects (red squares experimental data) 
the parameters are A = 0.62, a = 1.60, a' = 2.71, which lead 
to q = 1.37 and asymmetry a = —0.55. For negative defects 
(blue circles experimental data) the parameters are A = 0.58, 
a = 1.67, a' = 1.38, which lead to g = 1.49 and a — 0.14. 
The inset shows the peak region for positive defects. 


figure, these asymmetries can be fitted with the proba¬ 
bility densities introduced previously. 

In Fig. 6 we show a set of experimental data (e = 0.08) 
0, corresponding to the probability density p{vy) of 

the dimensionless velocity Vy = Vy/{Vy) — {vy)'^ of 

positive and negative defects (see also Fig. 4(a) in 
[llj). We find that these data can be very well fitted 
with the distribution Eq. (l30l) under the associations 

X ^ Vy - Aj v(5 = {Vyj ^{v'^) - {Vy)'^) - A/ -/d. Henco, 

p{vy) = aP{aVy - A/y/P), where a = sj (A^) - (A)2, 
follows from Eqs. ([28|) and ((29ll . Furthermore, we in¬ 
troduced an extra dimensionless parameter A that only 
shifts the complete distribution. Due to the previous 
rescaling, pivy) does not depend on the parameter /3. 
The parameters a, a', and A were determinate such that 
the global error is minimized. Even when the asymme¬ 
try is appreciable, the maximum of the distribution [Eq. 
(221 is around the origin. In fact, the influence of the 
shift introduced by A is only appreciable around the ori¬ 
gin (note that in both cases |A| < 1). For both positive 
and negative defects a very well fitting is obtained. We 
also checked that a similar fitting is obtained for higher 
values of the dimensionless temperature, e = 0.17, where 
the distribution asymmetry is smaller than that shown 
in Fig. 6 (see Fig. 4(b) in [ill). Hence, we conclude that 
dynamics of the defects velocity may be though as be¬ 
ing governed by two cascade mechanisms with different 
statistical properties, such as that defined by Eq. (1221) . 


T I ' I ' I ' I ' I ^ 
■ positive : 
> • negative 


-4 -2 0 2 4 i* 

.*! \ 

^ ft : 4 
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While a rigorous derivation of this interpretation is not 
developed here, the quality of the fitting gives a consis¬ 
tent support to the proposed theoretical frame. 

VI. CONCLUSIONS 

The present approach relies on expressing the variable 
of interest, associated to a given complex system, as a 
function of two independent Gamma random variables. 
These variables represent intrinsic fluctuations that drive 
the system. In addition, the complexity of the dynam¬ 
ics is represented by a random-system-response that also 
depends, in a non-linear way, on the fluctuations. Writ¬ 
ing the system response in terms of a generalized mean 
value, Eq. m, we showed that the arithmetic and geo¬ 
metric cases allow us to introduce a class of asymmetric 
g-Gaussian distributions. In the symmetric case, for any 
value of the complexity parameter g, they recover den¬ 
sities that follow from Tsallis entropy maximization. A 
similar approach applies for g-exponential distributions, 
which become defined in terms of a distortion parameter. 
We also showed that the complete family of obtained dis¬ 
tributions can be related via a change of variables with 
a Beta distribution. A g-triplet was derived for the sym¬ 
metric and unmodihed distributions. 

On one side, these results define an alternative numeri¬ 


cal tool for random number generation obeying the previ¬ 
ous statistical behaviors. On the other hand, the present 
approach may provide an alternative and very simple ba¬ 
sis for understanding statistical behaviors in complex dy¬ 
namics. Of special interest is the possibility of relating 
any asymmetry in the probability distributions with dif¬ 
ferent underlying cascade mechanisms. We have shown 
that in fact asymmetric Poissonian g-Gaussian densities 
(1 < g < 3) provide a very well fitting to the statisti¬ 
cal distribution of log-return signals in financial markets 
(Fig. 5) as well as the probability distribution of the ve¬ 
locity of moving defects in inclined layer convection m 
(Fig. 6). Therefore, the derivation of the present ap¬ 
proach from deeper microscopic or mesoscopic descrip¬ 
tions is an issue that with certainty deserves extra anal¬ 
ysis. The possibility of recovering asymmetric distribu¬ 
tions or the Beta statistics from non-extensive thermo¬ 
dynamics also remains as an open problem. 
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